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ABSTRACT 

Microwave emission of solar flares is formed primarily by incoherent gyrosynchrotron radiation 
generated by accelerated electrons in coronal magnetic loops. The resulting emission depends 
on many factors, including pitch-angle distribution of the emitting electrons and the source 
geometry. In this work, we perform systematic simulations of solar microwave emission using 
recently developed tools (GS Simulator and fast gyrosynchrotron codes) capable of simulating 
maps of radio brightness and polarization as well as spatially resolved emission spectra. A 3D 
model of a symmetric dipole magnetic loop is used. We compare the emission from isotropic 
and anisotropic (of loss-cone type) electron distributions. We also investigate effects caused by 
inhomogeneous distribution of the emitting particles along the loop. It is found that effect of 
the adopted moderate electron anisotropy is the most pronounced near the footpoints and it also 
depends strongly on the loop orientation. Concentration of the emitting particles at the loop top 
results in a corresponding spatial shift of the radio brightness peak, thus reducing effects of the 
anisotropy. The high-frequency (> 50 GHz) emission spectral index is specified mainly by the 
energy spectrum of the emitting electrons; however, at intermediate frequencies (around 10-20 
GHz), the spectrum shape is strongly dependent on the electron anisotropy, spatial distribution, 
and magnetic field nonuniformity. The implications of the obtained results for the diagnostics of 
the energetic electrons in solar flares are discussed. 

Subject headings: radiation mechanisms: non-thermal — Sun: corona — Sun: flares — Sun: radio radiation 



1. Introduction 

Microwave emission produced during solar 
flares is known to contain highly important infor- 
mation about fast electron acceleration and trans- 
port, coronal magnetic field and thermal plasma 
(jBastian et~aHll998tlFleishman et~al]|2009l) . How- 
ever, this potential of the microwave emission has 
not yet been converted to a routine diagnostics 
because of two main reasons. The first of them 
is absence of well-calibrated radio imaging spec- 



troscopy data with n eeded s patia l, temporal, and 
spect ral resolutions (|Garvl l2003t iGarv fc Keller! 
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2004). Currently, the situation has started to 
change as a number of solar radio instruments 
(e.g., Owens Valley Solar Array, OVSA, and 
Siberian Solar Radio Telescope, SSRT) experience 
significant upgrade (Upgraded SSRT, USSRT) and 
expansion (Expanded OVSA, EOVSA), while even 
more powerful solar radio facilities are planned to 
be built in near future, and a general-purpose ra- 
dio facility, EVLA, will soon be operational in 
the solar observing mode. This implies that the 
required radio data will soon become available. 
However, even with them there is a second reason, 
which is apparent lack of realistic 3D modeling of 
the microwave emission from flares. This modeling 
is highly important because the gyrosynchrotron 
(GS) emission depends in a complicated nonlin- 
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ear way on many involved parameters and source 
geometry including spatial inhomogeneity and an- 
gular anisotropy. The realistic modeling has to 
establish a clear quantitative picture and solid de- 
tailed understanding of how the involved physics 
(i.e., source properties and parameter regimes) af- 
fects the emission produced, e.g., what changes in 
the emission can be expected from variation of a 
given parameter 



on the radiation in the coronal plasma including 
polarization modification, which might be present 
in the reality, is ignored. A reasonably quick 
computation of a given model is made possible 
by the use of recently develo ped fast GS codes 
( Fleishman fe Kuznetsov 2010l) that proved to be 
highly accurate for both isotropic and anisotropic 
electron distributions. We intentionally consider 



here analytical models for the magnetic field and 
Ava ilable 3D mod els of th e GS emission (IPreka-Papaden^egtr An^rjjp^in to fully control the input and 
1993 iKucera et al.l Il993t iBastian et alj Il998l rel iably interpret the outcome; numerical input in 

Tzatzakis et a l. 2008: Simoes fe Costal200alFleishman et fe rm of corresponding datacubes is under develop- 
bimoes stal20inF have established valu- 



2009; 



able examples of the flaring microwave emis- 
sion; however, they are neither numerous nor 
comprehensive and rely on isotropic pitch-angle 
distribution and uniform spatial distribution in 
most of the cases. The observations, however, 
suggest that fast electrons in flares have often 
aniso tropic and/or inhomogeneous distributions 
fe.g. jMelmkov et al.ll2002t iFleishman et al.|[200l 
FleishmanlbOOd lAltvntsev et al.ll2008l ). 



Analysis of the pitch-angle anisotropy effect 
has yielded controversial conclusions: although 
the anisotropy was found to have huge effect on 
the GS emission from a uniform (or spatial! 
resolved) source ( Fleishman fc Melnikovl 12003 



Fleishman et al.1 |2003), this effect can become 



much weaker when averaging over significant vol- 
ume with a nonuniform loop mag netic field comes 
into play (|Simoes fc Costal 120101 ) . This calls for 
a more systematic study of the GS emission 
from anisotropic electron distributions in 3D mag- 
netic loops. The situation is additionally compli- 
cated by the fact that the mentioned pitch-angle 
anisotropy implies, as a side effect, a spatial in- 
homogeneity of the electron distribution due to 
the fast electron accumulation at t he loop-top 
(|Lee et alJll994t iMelnikov et alJl2002h . 

As a first step towards addressing the whole 
problem of GS modeling in realistic 3D coro- 
nal magnetic loops, this paper presents a conve- 
nient modeling tool, "3D GS Simulator", which 
gives any user an ability to build an analytical 
dipole magnetic loop, select a desired viewing an- 
gle, populate the loop with thermal plasma and 
nonthermal electrons, and compute microwave im- 
ages and spatially resolved spectra. The vacuum- 
like propagation from the source to the observer 
is adopted in the tool, so any propagation effect 



ment and will be presented soon elsewhere. This 
paper considers the microwave (GS + free-free) 
emission produced by mildly anisotropic electron 
distributions from a symmetric dipole magnetic 
loop viewed at different angles and for different 
parameter combinations. We find that even in 
such non-extreme conditions the anisotropy makes 
noticeable fingerprints on the emission properties, 
which are discussed in detail below. 

2. Simulation tool and method 

As has been said, the diversity of the microwave 
flaring emission is huge, so neither a paper nor a 
paper series can, perhaps, offer a truly compre- 
hensive table of options fully covering all relevant 
parameter combinations. Thus, entirely different 
approach to address the whole problem of the GS 
modeling is called for: we need a simulation tool 
capable to smoothly change all the involved source 
parameters and quickly compute and return the 
datacubes describing the microwave emission pro- 
duced. Here, we present such a tool and give an 
example of using it for the microwave emission 
simulation. 

2.1. Dipole magnetic flux tube model 

The simulations were performed using the in- 
teractive IDL tool GS Simulator. This tool allows 
one to change the shape and orientation of the flar- 
ing loop, choose the parameters of the magnetic 
field, thermal plasma, and energetic electrons, and 
calculate the brightness maps of GS emission. In 
the model used, the magnetic field is produced by 
a dipole located below the solar surface. The loop 
is formed by a set of field lines such that at the 
loop top it has a circular cross section with a given 
radius; near the footpoints, the loop becomes nar- 



2 



(a) 



rower and non-circular due to conservation of the 
magnetic flux. As shown in Fig. [1] which pictures 
an actual implementation of such geometry in the 
GS Simulator tool, the user is allowed to freely 
rotate the dipole loop model in any direction, so 
that to obtain an arbitrary line of sight relative to 
the dipole's central plane. 

The magnetic model geometry and adjustable 
parameters are described in detail in Appendix. 
The source code and documentation of the GS 
Simulator tool are provided as an online supple- 
ment. 




Fig. 1. — The magnetic field model used in the 
simulations, as implemented in the GS Simulator 
tool. The dashed box inscribes the portion of the 
magnetic loop (visualized by a few bold held lines) 
situated above the solar surface, while the solid 
rectangle represents the top view of an inscribing 
box that is perpendicular to the observer's line of 
sight. The two panels show two different orienta- 
tions of the same model corresponding to a loop 
located near the center of the solar disk (a) and a 
loop located near the solar limb (b). 



2.2. Flaring loop model and simulation 
method 

To keep the range of options manageable within 
one paper, we restrict in the present study the 
flexibility of the model parameters in several ways. 
We consider only a simple case of symmetric dipole 
magnetic loops, i.e. with the dipole perpendic- 
ular to a local vertical. The adopted geometry 
(visualized via a few reference magnetic held lines 
demarcating the surface of the magnetic loop) for 
two different loop orientations is shown in Fig. [T] 
Dependence of the magnetic field strength at the 
loop axis on the distance from the loop top along 
the held line is shown in Fig. 



We assume that the density and temperature of 
the thermal plasma component within the loop are 
constant (since the plasma in flaring loops is often 
heated up to the temperatures of > 10 7 K, the cor- 
responding barometric scale heights far exceed the 
loop heights, so the density variations with height 
can be neglected). Parameters of the energetic 
electrons can either be constant or vary with the 
distance from the loop top (see below). By using 
the above assumptions, we construct a 3D model 
of the flaring loop (which is observed at a given 
direction) with all the source parameters depend- 
ing on the Cartesian coordinates (x, y, z), where x 
and y are the coordinates in the image plane, and 
z is the coordinate along the line of sight. The 
source volume is divided onto a number of vol- 
ume elements (voxels); each voxel is considered to 
be quasi-homogeneous. The radio brightness map 
(or the observed emission intensity as a function 
of 2D coordinates x and y) at a given frequency / 
is calculated by numerical integration of the radi- 
ation transfer equation 



&h,n.{f,x,V,z) 



dz 



Jl,r(/, x, y, z) 



along all selected lines of sight. In Eq. (JTJ, II 
and 2r are the spectral intensities of the left- and 
right-polarized emission components, respectively, 
j'l and jr. are the corresponding GS emissivities, 
and and % are the absorption coefficients. We 
use the "weak coupling" model, i.e., the left- and 
right-polarized emission components propagate in- 
dependently. Left-polarized emission corresponds 
to either ordinary or extraordinary magnetoionic 
mode, depending on the magnetic field direction; 
respectively, right-polarized emission corresponds 
to the opposite mode. The plasma emissivities 
jo,x and absorption coefficients xo,x for the or- 
dinary and extraordinary modes accounting for 
both GS and free-free contributions at each voxel 
are calculated using fast GS c odes developed by 
Fleishman fc Kuznetsov ( 2010t ). Outside the flar- 
ing loop, the emission propagates like in a vacuum. 

The loop orientation is described in general by 
three Euler angles; however, since rotation of the 
loop around z axis results simply in the same ro- 
tation of the brightness maps, variation of two an- 
gles only is considered. We adopt the loop to be 
located at the solar equator; in this case, the loop 
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Fig. 2. — Magnetic field strength (a), loss-cone boundary (b), and relative density of the energetic electrons 
(c) vs. coordinate along the loop. All the values correspond to the loop axis. 



orientation is characterized by the angle ip between 
the magnetic dipole and the equatorial plane and 
by the longitude A. 

The energetic electrons are described by the 
distribution function F in a factorized form: 
F(E,fi) = u(E)g(fi), where E is the electron 
kinetic energy, /i = cos a, and a is the electron 
pitch-angle (the angle between the particle veloc- 
ity and the local magnetic field vectors). The 
tool allows using a number of different model 
electron distribution functions, including ther- 
mal, power-law, broken power-law, a power-law 
matched to a maxwellian core (so-called, ther- 
mal/nontermal, TNT, distribution), kappa dis- 
tribution, etc. In particular, use of TNT distri- 
bution allows one to consistently take into ac- 
count both GS (nonthermal) and gyroresonant 
(GR, thermal) contributions to the emission and 
absorption. For our modeling, however, it has a 
disadvantage that two distributions with different 
anisotropics must be matched, which eventually 
increases the number of model parameters needed 
to be varied in the modeling. For this reason, we 
assume that the electrons have a power-law en- 
ergy spectrum u(E) ~ E~ in the energy range 



EL 



< E < E n 



and so neglect the GR contri- 



bution entirely. However, we make an assessment 
of the GR contribution by considering the TNT 
distribution in the on-line supplement album. Al- 
though the effect of GR absorption is modest for 
our adopted source model and parameters, it does 
have a noticea ble imprint on the spectrum and 
polar ization (sf. lPreka-Papadema fc Alissandrakis 
Il992l ) from the lower parts of the loop legs, see 
Figs. 29-30 of the on-line album. 

The pitch-angle distribution can be either 



isotropic or a loss-cone described by the model 
function 



g(p) 



1, 

exp 



(ImI-Mc) 2 



for \fi\ < 
for \fi\ > fi c 



(2) 

where fi c = cos a c , a c is the loss-cone boundary, 
and the parameter A/z determines the sharpness of 
this boundary. The loss-cone boundary is adopted 
to exactly follow the transverse adiabatic invariant 



sin 2 a c 



B 



(3) 



where B and B[ are the magnetic fields at a given 
point and at the loop footpoint, respectively. De- 
pendence of the loss-cone boundary on the coor- 
dinate along the loop is shown in Fig. [2Jd; this 
parameter equals 90° at the footpoint and de- 
creases with height, so that the distribution be- 
comes closer to the isotropic one. Since within the 
adopted model the anisotropy is only strong at 
and around the footpoints, while becomes much 
weaker across most of the magnetic loop, we call 
this "a moderate anisotropy" . 

We consider both the homogeneous spatial dis- 
tribution of the energetic electrons (when their 
number density n c is constant) and the case when 
energetic electrons are accumulated at the loop top 
( Melnikov et al.ll2002h . The inhomogeneous distri- 
bution of the particles along the loop is described 
by the following model function: 



exp[-e 2 (^-7r/2) 2 ] 



(4) 



where tp is the magnetic latitude (or the angle be- 
tween the dipole axis and the vector drawn from 
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the dipole center to a given point) and the param- 
eter e determines the inhomogeneity degree (e = 
corresponds to the homogeneous case). Density 
profiles of the energetic electrons along the loop 
for different values of e are shown in Fig. 

In the simulations, we use the following pa- 
rameters of the flaring loop: height H = 10" ~ 
7 270 km from the base of the corona, dipole 
depth below the base of the corona D — 6" ~ 
4 360 km (so that the distance between the foot- 
points A ~ 11.5" ~ 8 400 km), radius at the top 
i?t = 2" ~ 1 450 km, and magnetic field at the 
top B t = 75 G (that results in the footpoint mag- 
netic field of Bf ~ 800 G). Figures UK] correspond 
to these parameters. Two loop orientations are 
considered in the paper in some detail: tp = 60°, 
A = 20° (a loop near the center of the solar disk, 
Fig. CQi) and ip = 60°, A = 80° (a loop at the 
limb, Fig. [TJd); many more examples are given 
in the on-line album. The thermal plasma den- 
sity and temperature are uq = 10 10 cm -3 and 
To = 2 x 10 7 K, respectively. The energetic elec- 
trons have the power-law index S — 4, cutoff en- 
ergies E m i n — 100 keV and i? max = 10 MeV, and 
the loss-cone boundary width A/i = 0.2. Thus, 
in each loop orientation, the variable parameters 
are: type of the pitch-angle distribution (isotropic 
or loss-cone), the number density of the energetic 
electrons n c , and the inhomogeneity parameter e. 

3. Simulation results 

3.1. Effect of the electron anisotropy 

Firstly, we consider effect of electron anisotropy 
on the GS emission. Figure [3] shows the brightness 
maps at four frequencies for the loop located near 
the center of the solar disk. The displayed Stokes 
parameters are: I = Jr+Jl and V = /r — Jl. The 
isotropic (top row) and loss-cone (bottom row) 
pitch-angle distributions are considered. In both 
cases, the number density of the energetic elec- 
trons n c is constant and equals 3 x 10 6 cm" 3 . Fig- 
ure[4]shows the spatially resolved emission spectra 
for the footpoint and loop top sources (obtained 
by summation over the circled pixels) as well as 
the total emission (obtained by summation over 
all pixels) . The figure also shows the degree of po- 
larization r\ and spectral index <5 r which are defined 



V= T ^ = -JJJ- (5) 

From the spatially resolved spectra, one can 
see that the influence of the electron anisotropy 
is most important at the footpoints. In the opti- 
cally thin frequency range, the emission from the 
loss-cone distribution is lower than that from the 
isotropic electrons by a factor of about 2-6. This 
reflects the fact that the GS radiation is emit- 
ted mainly in the direction of the electron veloc- 
ity. Near the footpoints, the electrons with the 
loss-cone distribution are concentrated around the 
pitch-angle of 90° whereas the magnetic field (in 
the adopted geometry) is nearly parallel to the line 
of sight. This is why the electrons with the loss- 
cone distribution produce only a relatively weak 
radiation flux towards the observer. In the op- 
tically thick frequency range, the emission inten- 
sities from the isotropic and anisotropic distribu- 
tions are almost the same, although the differences 
in the degree of polarization and spectral index 
can be visible. Near the loop top, the loss-cone 
boundary a c falls to about 20° and thus the loss- 
cone distribution does not noticeably differ from 
the isotropic one; as a result, these distributions 
produce almost identical emission. 

A complementary way of thinking of the emis- 
sion is via the radio images at various frequencies 
(Fig. U). At the low frequencies (3.75 GHz), the 
whole loop is seen on the map (although the foot- 
points are brighter than the top); the images for 
the isotropic and anisotropic distributions are very 
similar. At the higher frequencies, the emission is 
strongly concentrated at the footpoints. However, 
since the footpoint emission from the electrons 
with the loss-cone distribution is weaker than that 
from the isotropic distribution, the difference be- 
tween the footpoints and the loop top is smaller as 
well. Therefore the emission from the anisotropic 
electrons is more evenly distributed along the loop. 

Figures [5][6] show the brightness maps and emis- 
sion spectra for the loop located near the limb (all 
other parameters are the same as in Figures [3]|4|) . 
One can see that now the electrons with the loss- 
cone distribution produce stronger footpoint emis- 
sion in the optically thin frequency range than the 
isotropic distribution does (due to the same rea- 
sons as discussed above). In the optically thick 
frequency range, the considered characteristics of 
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Fig. 3. — Radio brightness maps for a loop located near the center of the solar disk, for the isotropic 
distribution (top row) and loss-cone distribution (bottom row); the north is up and the west is to the right. 
Concentration of the accelerated electrons is assumed to be constant along the loop. Solid lines are the 
intensity contours which are evenly distributed between zero and the maximum brightness temperature T m 
(the corresponding temperatures are given in each panel). Color shades (see online version of the journal) 
represent the circular polarization (Stokes V normalized by the absolute value of V peak); red and blue 
correspond to the right and left circular polarizations, respectively. 



the emission are very similar for both isotropic and 
anisotropic cases except the effects of the individ- 
ual cyclotron harmonics that are much more pro- 
nounced for the anisotropic distribution, which is 
discussed in Section [3T41 in greater detail. The loop 
top emissions from the isotropic and loss-cone elec- 
tron distributions are almost the same. 

In the brightness maps (Fig. [5]), the maximum 
of the emission at 3.75 GHz is located between 
the loop top and the footpoints; the images for 
the isotropic and anisotropic distributions are very 
similar. At the higher frequencies, the emission 
is concentrated at the footpoints. In contrast to 
the previous case (Fig. [3]), the footpoint sources 
are now more compact for the loss-cone electron 
distribution; this effect becomes more pronounced 
with the frequency increase. 

Changing the pitch-angle distribution from the 



isotropic one to the loss cone results in a shift 
of the spectral peak of the footpoint emission to- 
wards lower frequencies for the loop located near 
the disk center, and towards higher frequencies 
for the loop near the limb (which could be ob- 
servationally addressed via center-to-limb varia- 
tion analysis of the spatially resolved footpoint 
spectra). This spectral peak variation affects the 
emission spectral index around the peak. With an 
increasing frequency (in the optically thin range), 
the spectral indices of the emission from the 
isotropic and anisotropic electrons approach each 
other and gradually become the same. We have 
confirmed that at / — > oo, the emission spectral 
indices asymptotically approach the ultrarelativis- 
tic limit S te i — (5 — l)/2 (provided that the high- 
energ y cutoff -E max — > oo. lGinzburg fc SvrovatskiH 
I965L and neglecting the free-free contribution) 
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Frequency, GHz Frequency, GHz Frequency, GHz Frequency, GHz 

Fig. 4. — Emission intensity, degree of polarization, and spectral index vs. frequency for the loop shown 
in Fig. [3] The columns (from left to right) correspond to: northern footpoint source, southern footpoint 
source, loop-top source, and the emission from the entire loop (spatially unresolved). Solid lines: isotropic 
distribution; dashed lines: loss-cone distribution. The regions taken to calculate the spatially resolved spectra 
are indicated in Fig. [3] by thick dashed circles. 



for both the isotropic and anisotropic electron 
distributions. However, in the frequency range 
of ~ 10 — 100 GHz, the emission spectral index 
is strongly affected by the magnetic field inho- 
mogeneity in the source; as a result, the spectral 
index varies with frequency and can be noticeably 
different at the different footpoints and the loop 
top. At the frequencies > 100 GHz, the spec- 
trum is affected by the high-energy cutoff of the 
electron distribution (note that Figs. [3j6] are cal- 
culated for E max = 10 MeV) which makes the 
spectrum steeper (the spectral index increases 
with frequency); this effect is clearly visible for 
the loop-top emission where the magnetic field in- 
homogeneity is the lowest. On the other hand, at 
these and higher frequencies the free-free contri- 
bution can make the spectrum flatter (see below 



examples of such flattening). 

The loss-cone anisotropy of the fast electrons 
has opposite effect on the GS emission for the 
loops located near the solar disk center and the 
limb. Thus, for some intermediate longitude, the 
anisotropy effect should be minimal. We have 
found that this occurs at the longitude of A ~ 60°. 

3.2. Effect of the inhomogeneous electron 
distribution 

Figures [7]fTUl show the brightness maps and spa- 
tially resolved emission spectra for the models 
with the inhomogeneous distribution of the ener- 
getic electrons along the flaring loop (the particles 
are accumulated at the loop top). The number 
density of the energetic electrons at the loop top 
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Fig. 5. — The same as in Fig. [31 for a loop located near the solar limb. 



is taken to be n t = 2.8 x 10 s cm~ 3 . The density 
profiles of the energetic electrons are described by 
Equation (|4]) with e = 4, 6, and 8, which corre- 
sponds to the ratios of the footpoint and loop-top 
number densities of n,{/n t = 1.1 x 1CP 2 , 3.6 x 10 -5 , 
and 1.3 x 10~ 8 , respectively (for e = 4, we obtain 
the footpoint number density of rif = 3x 10 6 cm~ 3 , 
like in the previous Section). In all cases, the loss- 
cone distribution of the energetic electrons over 
the pitch-angle is used. 

Figures [7j|8] correspond to the loop located near 
the center of the solar disk. Since the number 
density of the energetic electrons at the loop top 
is assumed to be constant, the loop-top emission 
is almost independent on the parameter e. In the 
footpoints, the increase of e results in a decrease 
of the number density of the energetic electrons 
and, consequently, in a decrease of the emission 
intensity in the optically thin frequency range; the 
spectral peak shifts towards lower frequencies. In 
the brightness maps, the intensity maximum is al- 
ways located near the loop top. For e = 4, at the 
frequency of 3.75 GHz, the whole loop is seen with 
relatively weak brightness variations along it. At 



higher frequencies and/or with increasing e, the 
visible emission source becomes more compact. 

Figures IMTOl show the brightness maps and 
emission spectra for the loop located near the limb. 
Like in the previous case, the loop-top emission is 
independent on the parameter e while the foot- 
point emission intensity (in the optically thin fre- 
quency range) decreases with increasing e. Dis- 
tributions of the radio brightness along the loop 
can be qualitatively different for the different fre- 
quencies and inhomogeneity models: for e = 4, 
the emission maximum at 3.75 GHz is located at 
the loop top and gradually shifts towards the foot- 
points with increasing frequency. For e = 6 and 
8, the intensity maximum is always located at the 
loop top, but the visible emission source becomes 
more compact with increasing frequency. Thus at 
high frequencies (> 15 GHz), a brightness map 
will be dominated by either two footpoint sources 
(for the loops with a moderate accumulation of 
the particles at the top, up to n%/nf ~ 100) or 
one loop-top source (for the loops with higher in- 
homogeneity). We note an interesting behavior of 
the polarization map: for the loop-like optically 
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Fig. 6. — The same as in Fig. HI for the loop shown in Fig. [5] (located near the solar limb). The regions 
taken to calculate the spatially resolved spectra are indicated in Fig. [5] by thick dashed circles. 



thin sources (at 17 or 34 GHz) in Figure |H] most 
of the image area is dominated by the RCP sense, 
which is indicative of the same direction of the 
line-of-sight magnetic field component throughout 
most of the loop area. Similar unipolar polariza- 
tion patterns are highly typical for the NoRH data 
on the limb flares. If the concentration of the 
electrons at the loop top develo ps at the course 
of flare due to electron transport (ILee et al.lll994 ; 



Melnikov et all 120021; iFleishmanl 120061 ). this pro- 



cess of the electron inhomogeneity build-up can 
be roughly described in terms of the e increase, 
which implies the apparent brightness peak mo- 
tion (at the optically thin frequencies) from foot- 
points to the looptop in ag reement with observa- 



tions of some limb fl ares (jTzatzakis et al 
Reznikova et al. 20091 ). 



2008; 



In all the models, the emission spectral index at 
the loop top behaves like that for a homogeneous 
source: the index is nearly constant (c5 r ~ 1.7—1.8) 



in the frequency range of 20-50 GHz (for both loop 
orientations), and increases at higher frequencies 
due to the high-energy cutoff of the electron spec- 
trum. At the footpoints, an increase of the param- 
eter e results in a decrease of the spectral peak 
frequency. With an increasing frequency (in the 
optically thin range), the spectral indices for all 
values of e become the same (if the free-free emis- 
sion is negligible, see below). However, in con- 
trast to the loop top, the spectral index of the 
footpoint emission is also affected by the source 
inhomogeneity. For the loop located near the disk 
center, we can notice an asymmetry of the emis- 
sion spectra produced at the different footpoints, 
which is caused by a difference of the viewing an- 
gles. For the loop located near the limb, the high- 
frequency spectral index of the footpoint emission 
for the model with e = 8 differs noticeably from 
the indices for e = 4 and 6. This is caused by 
the contribution of the free-free emission of the 
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Fig. 7. — Radio brightness maps for a loop located near the center of the solar disk for the different 
distributions of the accelerated electrons along the loop. The accelerated electrons are assumed to have the 
loss-cone distribution. 



thermal electrons which becomes dominant under 
these conditions: at large frequencies the spectral 
index of the free-free emission goes to zero. 

Changing the pitch-angle distribution of the 
energetic electrons affects the spatially resolved 
emission spectra in the same way as discussed 
in the previous Section: in comparison to the 
isotropic distribution, the loss-cone one provides 
stronger footpoint emission for a loop near the 
limb and weaker footpoint emission for the loop 
near the disk center, while the loop-top emission 
remains almost unchanged. Since in the consid- 



ered inhomogeneous models the loop-top emission 
source often dominates the images, the effect of 
the electron anisotropy on the brightness maps is 
only moderate. 

3.3. Spatially integrated spectra 

In this Section, we consider the total (spatially 
integrated) emission from the flaring loop. Al- 
though in the view of the currently operational 
and soon to be available radio instruments, consid- 
ering the total radiation may look old-fashioned, 
we feel that this still makes sense for the follow- 
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Frequency, GHz Frequency, GHz Frequency, GHz Frequency, GHz 

Fig. 8. — Emission intensity, degree of polarization, and spectral index vs. frequency for the loop models 
shown in Fig. [JJ Solid lines: e — 4; dashed lines: e — 6; dotted lines: e = 8. The regions taken to calculate 
the spatially resolved spectra are indicated in Fig. |3]by thick dashed circles. 



ing reasons. Most of the historically accumulated 
databases and corresponding statistical studies are 
done based on the total power observation s (e.g., 



Guidice fc Castellil 1197.4 iNita et all 120041 ) . The 



total power data are simpler manageable as they 
can be easily visualized by dynamic spectra and 
characterized by only a few simple numbers, such 
as spectral indic es, rise and decay times, peak flux 



and frequency (jNita et al.l 120041 ) . In particular, 



the corresponding high-frequency spectral index 
is widely used t o evaluate the fast ele ctron energy 
spectral index ( Dulk fc Ma rsh 1982 ). From this 
perspe ctive, a statement made by lSimoes fc Costal 
(|2010h that such a spectral index is a good measure 
of the electron spectral index even for anisotropic 
electron distributions, if confirmed, could be of a 
great practical value. 

The parameters characterizing the total emis- 
sion are shown in the right columns in Figs. [4j [6j [Sj 



and[T0] Considering the right top panel in Figure|4] 
as a vivid example, one can easily isolate three dis- 
tinct regions of the spectra — low-frequency part 
(region I), middle- frequency part (region II), and 
high-frequency part (region III). Visual compari- 
son of these total power spectra with the spatially 
resolved spectra from the footpoints and looptop 
in Figure 0] suggests that the low-frequency part 
is formed primarily at the looptop region with low 
magnetic field, the high-frequency part in the foot- 
points where the magnetic field is large, while the 
middle-frequency part by the entire loop and so 
related to the magnetic field non-uniformity. 

The low-frequency part is known to be formed 
by the effect of the GS optical thickness and/or 
the Razin-effect (suppression of th e GS emission 
in a d ense background plasma, e.g.. iBastian et all 
1998) possibly accompanied by th e free-free ab- 
sorption in the dense plasma f Ba stian et al 
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X, arcsec X, arcsec X, arcsec X, arcsec 

Fig. 9. — The same as in Fig. [7J for a loop located near the solar limb. 



20071 ); the slope of the spectrum can here be 
quantified by the index of —2 or less, see the 
right bottom panel in Figure 2) The high- 
frequency part is mainly determined by the dis- 
tribution o f fast electrons incl uding the energy 
spectrum ( Bastian et al.l 1998 ) and pitch -angle 
anisotropy ([Fleishman fc Melnikovl l2003al) ; note 
the anisotropy-rclatcd difference between the solid 
(isotropic) and dashed (loss-cone) curves in this 
panel. 

The middle-frequency part is clearly seen in 
Fig. S] as it is almost flat (the spectral index is 
around zero). Solar microwave bursts with flat 
spectra have been observed for decades. For ex- 



ample, iRamatv fc Petrosianl (|1972n proposed that 
the free-free absorpt ion of GS emissi on can form 
such spectra, while iLee et al.l ([19941 ) recognized 
that the electrons trapped in a large dipole mag- 
netic loop can produce the flat radio spectra due 
to the source non-uniformity in certain parame- 
ter regi mes. Our modeli ng confirms the finding 
made by ILee et al. I ([19941) . In addition, we have 
found that this middle-frequency part is not al- 
ways flat but can have either negative or positive 
spectral index, see, e.g., Fig. [51 which can be mis- 
interpreted as either region I or III in observations 
with a limited spectral coverage . In fact, the ob - 
servations (see, e.g., Fig. 11 in iNita et al] |2004[ ) 
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Frequency, GHz Frequency, GHz Frequency, GHz Frequency, GHz 

Fig. 10. — The same as in Fig. HI for the loop models shown in Fig. |9] (the loop is located near the solar 
limb). The regions taken to calculate the spatially resolved spectra are indicated in Fig. [5] by thick dashed 
circles. 



reveal that the histograms of both low- and high- 
frequency spectral indices extend to zero imply- 
ing that both low- and high-frequency spectra can 
be much flatter than those determined by optical 
thickness effect or electron energy index, respec- 
tively. For practical application, this means that 
having the spectrum falling with the frequency 
does not guarantee that its slope is formed by ei- 
ther energetic or angular properties of the electron 
distribution function but can instead be related to 
the source non-uniformity. 

One of the parameters controlling the shape of 
the middle-frequency part is the viewing angle: 
Figure [5] presents the case when this part grows 
with frequency, which overall broadens the spec- 
trum peak. Not surprisingly, non-uniform spatial 
distribution also significantly affect this part of the 
spectrum. In fact, with an inhomogeneous non- 
thermal electron spatial distribution (with their 



concentration at the looptop) , the radio spectrum 
begins to resemble emission from a roughly uni- 
form source, see Figures [5] and [TO] The reason 
for this to happen for an inhomogeneous source 
is very simple: with the adopted inhomogeneous 
electron distribution most of them reside at the 
looptop, where the spatial variation of the mag- 
netic field and the viewing angle are minor, so we 
have a situation similar to a uniform source. 

Finally, let us consider to what extent we could 
rely on the high-frequency spectral index in eval- 
uating the fast electron energy spectrum or pitch- 
angle anisotropy. The obtained spectral behavior 
of the local spectral index can easily be tracked 
in the figures, so we do not describe it here in 
any detail. Instead, in order to quantify it with 
measures relevant to available observations, we 
form a few distinct spectral indices, namely, an in- 
dex Sis at 15 GHz relevant to the OVSA spectral 
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Table 1 

Spectral parameters of the total (spatially unresolved) emission. 





c 


/peak, GHz 


Sis 


#17-34 


#50 


Loop location 


: near the disk center 






Isotropic 





6.03 


2.40 


2.46 


2.05 


Loss cone 





4.37 


2.98 


2.21 


1.86 


Isotropic 


4 


7.59 


1.75 


1.82 


1.79 


Loss cone 


4 


7.59 


1.76 


1.77 


1.77 


Loss cone 


6 


6.92 


1.66 


1.72 


1.76 


Loss cone 


8 


6.92 


1.64 


1.71 


1.75 


Loop location 


: near the limb 








Isotropic 





12.02 


0.57 


2.10 


2.14 


Loss cone 





13.80 


0.13 


1.67 


2.16 


Isotropic 


4 


10.47 


0.94 


1.87 


1.92 


Loss cone 


4 


10.47 


0.79 


1.72 


1.95 


Loss cone 


G 


7.24 


1.77 


1.78 


1.78 


Loss cone 


8 


6.61 


1.66 


1.71 


1.76 



Note. — The values 615 and ($50 are the exact spec- 
tral indices at 15 and 50 GHz, respectively; #17-34 is 
an approximate spectral index calculated using two 
points at 17 and 34 GHz; for comparison # re i = 1.5 
while #d-m = 2.38. For the models with e = and 
e ^ 0, the number densities of the energetic electrons 
at the loop top are different (see text). 



range (|Nita et al.ll2004l )- the "Nobeyama" index 
#17-34 calculated with two well separated frequen - 



Reznikova et al. 2009), 



cies, 17 and 34 GHz (e.g. 
and higher-frequency spectral index #50 taken at 
50 GHz relevant to a number of earlier separated- 
frequ ency observations (e.g., Melnikov k, Magunl 
19981 ); the results are gathered in Table [T] 

One can see, that for the same energy elec- 
tron distribution, # — 4, the measured spectral 
indices are highly different depending on the pitch- 
angle anisotropy, the viewing angle, spatial inho- 
mogeneity of the fast electron distribution, and 
the instrument used to measure the spectral in- 
dex. In fact, the numbers given in Table [T] are 
bounded between 0.13 and 2.98, although as has 
already been said, the ultrarelativistic approxi- 
mation predicts the emission spectral index of 



#rci = (# — l)/2, or 8 Ie ] = 1.5 for the adopted 
value of 8 = 4, while a widely used ap proxima- 
tion proposed by iDulk fc Marsh (1982) predicts 
the value of # D -m ^ 0.908 - 1.22, or # D _ M ^ 2.38 
for 8 — 4; none of this approximate values is fa- 
vored by Table [TJ Stated another way, the radio 
indices from 0.13 to 2.98 would imply the electron 
energy index range from 1.26 to 6.96 if the rela- 
tivistic asymptote is used and from 1.50 to 4.67 if 
the Dulk-Marsh approximation is applied. Again, 
none of the values is truly informative in terms 
of recovering the original electron energy index as 
there is no unique way of linking the measured ra- 
dio spectral index with the original electron energy 
spectral index. This conclusion holds (with nar- 
rower scatter of the recovered values, 4.34 — 5.92 
and 3.21 — 4.09, respectively) if we consider only 
higher- frequency indices #17-34 and 850. On the 
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other hand, the OVSA index S15 is more sensi- 
tive to the particle anisotropy, which, thus, can be 
studie d through a forward fit of the radio spectra 
(l2009h . 



like in iFleishman et al 



3.4. Harmonic structure 

The GS emission from a homogeneous source 
can demonstrate an oscillatory spectral structure 
in the low- frequency range (///b < 10), when 
the emission intensity increases at a narrow spec- 
tral region near the harmonics (small integer mul- 
tiples) of the electron cyclotron frequency (see, 



e.g., the figures i n the ar ticles of|Ramat 
Benka fe Holmanl 11992 ; IFleishman fc Mclnikov 



3Y (see, 
vlll969t 



2003allbl; IFleishman fc Kuznetsovl l2010h . In an 
inhomogeneous source, however, this harmonic 
structure can be hidden because of natural smooth- 
ing: the resonance giving rise to a gyroharmon- 
ics at a given location will vary with frequency 
due to the spatially dependent resonant condi- 
tion in the spatially non-uniform magnetic field. 
Thus, even if a spectrum from a single pixel con- 
tains harmonics they often disappear after inte- 
gration over even a relatively small part of the 
source. This is why no harmonic structure is 
present in either footpoint or looptop spectra 
in Figures 0J El and [10] other than a num- 
ber of extremely large peaks in Figures [6] and 
ITTJ1 for the loss-cone case provided that the con- 
ditions for the electr on cyclotron maser (ECM) 
insta bi lity (see, e.g. 
19791: iHolman etal 



Stepanovl Il978l IWu fc Lee 
1980 ; iHewitt et all Il98a 



Dulk fc Marshlll982t ISharma fc Vlahosjll984l: IWu 



19851: lAschwandenl Il990l: IFleishman fc Melnikov 



199^ IFleishman fc ArznerlbOOOHLaBelle fc Treumann 



200S 



Treumann 20061: Kuznetsovl 12011 ) are lo 



cally fulfilled. We do not consider the coherent 
ECM emission here concentrating on the harmonic 
structure of the incoherent GS emission. 

In the simplest case of a uniform GS source, 
the harmonic structure is more prominently pro- 
nounced for the source viewed at a quasitransverse 
direction relative to the source magnetic field. 
Note, that for the magnetic model adopted here, 
the harmonic structure is only expected from the 
footpoint and leg regions, while not from the loop- 
top where the magnetic field is too weak for the 
gyroharmonics to be produced at the considered 
parameter regime; thus, the limb location of the 
loop is most favorable to produce distinct gyro- 



harmonics, which is confirmed by comparing Fig- 
ures unii 

These figures show the brightness maps at low 
frequencies for two loops locations — at the disk 
center and limb. As expected, there is no appar- 
ent harmonic structure in the images seen at the 
disk center, Fig. [TT] In contrast, for the loop lo- 
cated near the limb, the images contain a number 
of bright stripes highlighting the isolines of the 
magnetic field strength corresponding to the gy- 
roresonance conditions at a given frequency. At 2 
GHz, for example, these stripes are clearly visible 
from the footpoints up to about half of the loop 
height. If one were gradually increasing the fre- 
quency, for which the image has been computed, 
each stripe would demonstrate an apparent down- 
motion because same harmonic number requires 
proportionally larger magnetic field for a higher 
frequency. Accordingly, at higher frequency im- 
ages, the stripes are shifted toward the footpoints, 
and their contrast (or amplitude of the intensity 
oscillations) decreases. 

Apparently, having these stripes detected in 
real observations would offer a highly efficient way 
of model-independent quantitative measurement 
of the coronal magnetic field (note that averag- 
ing along the line of sight plays here only a minor 
role since the GS optical depth is large at the gyro 
harmonics so the layer contributing to the emis- 
sion is narrow implying almost uniform magnetic 
field along the essential part of the integration). 
In case of rather compact loop considered here 
this task requires a radio array with the spatial 
resolution about 1" at the low frequencies (where 
scattering of radio waves by coronal density in- 
homogeneities can additionally fuze the images, 
Bastianl (|1994l) . however), which is comparable 
with the anticipated capability of the Frequ ency 
Agile Solar Radiotelescope (TASR. lGarvll2003l ): al- 
though for the instruments currently under devel- 
opment (e.g., USSRT and EOVSA) detecting the 
harmonic stripes can only be expected from the 
biggest flaring loops and/or for the cases with a 
more uniform magnetic field. Such fa vorable cases, 
although untypical, are not unlikely: Staehli et al. 
(|1987n reported a total power data on a microwave 
burst obtained with a few single frequencies that 
demonstrated prominent harmonic-like enhance- 
ments at two frequencies. Furthermore, the OVSA 
total power data reveal clear horizontal stripes at 
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Fig. 11. — Radio brightness maps (emission intensity) for a loop located near the center of the solar disk. 
Brighter areas correspond to higher intensity. Number density of the accelerated electrons is constant and 
their pitch-angle distribution is of the loss-cone type. 
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Fig. 12. — The same as in Fig. [IT] for a loop located near the solar limb. 



the dynamic s pectra of some events, e.g. 
June 15 flare (jFleishman et al. 2008); see 



2003 
also 



(|Benka k Holmanlll 992). which might indicate the 



gyroharmonic contribution; however, such an in- 
terpretation remains ambiguous unless confirmed 
by direct imaging data. 

4. Discussion 

We have considered microwave emission pro- 
duced by moderately anisotropic electron distri- 
butions populating (uniformly or non-uniformly) a 
non-uniform symmetric dipole magnetic loop. We 
emphasize that for the adopted loop geometry, the 
mirror ratio is more than 10, so in most of the loop 
volume (except the close vicinity of the footpoints) 
the loss-cone angle defined by the adiabatic in- 



variant (Figure [5Jd) is small and, accordingly, the 
angular distribution of the fast electrons is close 
to the isotropic one, which we here convention- 
ally call a "moderate anisotropy". But even this 
moderate anisotropy affects noticeably both im- 
ages and spectral characteristics of the emission. 

In particular, for the uniform electron distri- 
bution along the loop the pitch angle anisotropy 
enhances the optically thin emission from the foot- 
point vicinity for the loop located at the limb, 
while suppresses it for the loop located at the disk 
center. In the case of non-uniform electron dis- 
tributions due to accumulation at the loop top, 
the relative contribution of the electrons in the 
vicinity of the footpoints decreases and, so, the 
effect of the anisotropy becomes less pronounced. 
Nevertheless, in most of the cases (except looptop 
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emission, where the electron angular distribution 
is almost isotropic) the effect of anisotropy on the 
images, spectrum, polarization, and spectral in- 
dex can easily be recognized. We note that for 
larger thermal density than adopted for our re- 
stricted modeling, w hich is indeed often the case 
( Bastian et al.l[2007h . the looptop emission will be 
strongly suppressed by the Razin effect (which is 
specified by the n /B ratio and becomes especially 
strong at / < 2Quq/B), so the emission from the 
footpoints and legs of the loop, where the electron 
anisotropy is stronger than at the looptop, will 
dominate the emission; thus the anisotropy effect 
will be even stronger than for the cases discussed 
above, see on-line supplement Album. 

An interesting feature of the low-frequency im- 
ages is the presence of the "gyro-stripes" , which 
are indicative of distinct gyroharmonics produced 
at certain heights changing with frequency (since 
a few small integer multiples of the local gy- 
rofrequency produce strong enhancement of the 
flux density). Having radio imaging spectroscopy 
data with the spatial resolution sufficient to re- 
solve such gyro-stripes would offer a nice model- 
independent diagnostics of the coronal magnetic 
field during flares. In the near future the Ex- 
panded OVSA will be capable of such measure- 
ments at least for relatively large flaring sources. 

An important issue is anisotropy /non- uniformi- 
ty effect on the total power radio spectrum. Since 
the total power spectrum is the result of emis- 
sion integration over the entire source area, the 
result of this integration depends on how the dis- 
tinct spatially resolved contributions are weighted, 
which, in turn, depends on the anisotropy, inho- 
mogeneity, and the viewing angle. We found that 
all the spectrum peak, peak flux, polarization, 
and the spectral index depend on both anisotropy 
and electron distribution inhomogeneity, although 
the effects are counter-directed. In particular, the 
radio spectral index changes noticeably with fre- 
quency and its behavior depends noticeably on the 
pitch-angle anisotropy, e.g., one can see that the 
maximal value of the total power spectral index 
(Figure |H right bottom) for the anisotropic case 
exceeds by roughly 1 that for the isotropic case. 

We compared the spectral indices with asymp- 
totic relativistic (synchrotron) and approximate 
Dulk-Marsh values as the observed radio spectral 
indices are often used to evaluate the fast elec- 



tron energy spectral indices. We found that nei- 
ther synchrotron nor the Dulk-Marsh value is a 
good approximation to the true spectral index: al- 
though the synchrotron index does represent the 
true index at high frequencies (higher than con- 
sidered in our plots), the Dulk-Marsh index is 
less meaningful even though it sometimes coin- 
cides with the true one at a given single frequency. 
Ther efore, we do not confir m the conclusion made 
by (|Sim5es k Costal l2010f) that the Dulk-Marsh 
index is a quantitatively good approximation to 
the spectral index of the spatially integrated spec- 
trum, which might be an artifact specific to their 
adopted source model. 

5. Conclusion 

We have introduced a flexible simulation tool 
capable of fast computing 3D models of the mi- 
crowave emission utilizing recently developed fast 
GS codes and presented an example of its use 
for the microwave emission modeling. Consid- 
ering a highly restricted parameter space (sym- 
metric dipole magnetic loop, moderate anisotropy, 
weak or no Razin-effect) we have analyzed im- 
ages, spectra, and polarization of the model ra- 
dio emission in the view of the available and fu- 
ture (more complete) radio data. In particular, we 
note that the high-frequency spectral index does 
not have a unique value for a given energy spec- 
tral index of radiating electrons; instead, it notice- 
ably varies with the frequency, the viewing angle, 
the anisotropy, and the inhomogeneity and also 
different for various parts of the source and the 
entire source. This implies that the use of the 
radio spectral index for constraining the electron 
energy spectral index is not that straightforward; 
in place, the electron spectrum recovery must rely 
on the forward fitting of the entire radio datacube 
(ideally, includin g polarization) like in t he exam- 
ple presented bv lFleishman et ah ( 20091 ). 
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Appendix: Magnetic Model 

Similar with other models used in the past 



(e.g. lAntiochos fc Sturrockl 119761 ) . in our simula- 
tion tool the magnetic field is produced by mag- 
netic dipole with a moment /x that makes an ad- 
justable angle n/2 — po with its corresponding so- 
lar radius, and located below the solar surface at 
a depth D. The flaring loop is constructed around 
a central field line that is chosen to lie in the plane 
defined by the magnetic dipole vector and the lo- 
cal vertical, referred from now on as the central 
plane of the loop. As pictured in Figure H31 in the 
dipole's cartesian system of coordinates C}; 
which is defined as having the axis 0£ oriented 
along the magnetic dipole fi, and the axis Or] in 
the central plane, the central magnetic field line 
may be parameterized as 



£ c = H sin 2 ip cos (p 
- H sin 3 pi, 



(6) 



where ip € [0, 180] represents the angle between 
the position vector r c = {£, r], 0} and the dipole's 
axis, i.e. cos<^ = £,/r c , and H is the height of 
the central line, i.e the maximum distance between 
the central line and the dipole's axis, which corre- 
sponds to (f = n/2. In a solar coordinate system 
defined as having the axis Oy along the solar ra- 
dius and the axis Oz identical with the axis O^, 
the same central line is expressed as 



H sin 2 (ip) cos (ip + ipo) 
H sin 2 (ip) sin (p + po) . 



(7) 



Equations ((6|) and Q are a direct consequence of 
the fact that any magnetic dipole field line is de- 
scribed by the general polar equation r = H sin 2 ip 
and that the dipole's system of coordinates is ro- 
tated relative to the solar system of coordinates 
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by the angle tpo around the Oz axis, i.e. 

x = £ cos ipo — r] sin tp 

y = £ sin p> + rj cos tpo (8) 

Z = C 

The magnetic flux tube centered on the central 
field line is defined in terms of a circular cross sec- 
tion of radius po that is normal to the central field 
at cp = 7r/2, i.e the magnetic loop-top. Hence, 
the cartesian coordinates of an arbitrary magnetic 
field line, r = {£,rj, (}, that intersects the loop- 
top cross section at distance p from the central 
field line, may be expressed in terms of three con- 
venient free parameters {p, a, ip} as 



= \JH 2 + 2pHcosa + 
= (H + pcosa) sin 3 ip 

• 3 

= p sin asm (p, 



p 2 sin 2 ip cos tp 



(9) 



where, as shown in Figure 113b . a represents the 
fixed angle between the vector p and the dipole's 
vertical axis Or], and p represents the variable an- 
gle between the vector r and the dipole moment fi 
that is measured in the plane containing the field 
line, which makes the angle 9 = arctan(C/??) = 
p sin a/ (H + p cos a) with the central plane of the 
magnetic flux tube. 

For p = po, the set of equations ((9J defines the 
envelope of the magnetic flux tube, which has a 
circular cross section shape only for tp = n/2, 
i.e. at the loop-top, while the shape of its cross 
section is continuously changing along the central 
field line. 

The condition for an arbitrary point r to be- 
long to a magnetic flux tube of height H and top 
circular cross section of radius po, i.e. the condi- 
tion for the magnetic field line passing trough the 
point {£, 7], C} to intersect loop-top circular cross 
section, may be expressed as 



H-r, 1 



e+v 2 



3/2' 



+ r i + 



e + if 



< 



which has to be combined with the condition 



x 2 + (y - D + Rsun? + z 2 > R 2 Sv 



(11) 



in order to determine if such flux tube point is also 
located above the solar surface. 



19 




Fig. 13. — The magnetic dipole moment model implemented in the simulation tool. 



Considering a normal cross section of the flux 
tube corresponding to an angle ip measured in the 
central plane of the flux tube, the parameter 



iff | cos yV 1 + 3 cos 2 <p + (12) 
COS if + v/l + 3cos 2 }, 



2 

-Lin 

V3 



may be used, instead of the angular coordinate ip, 
as a convenient flux tube longitudinal coordinate 
corresponding to the normal cross-section of inter- 
est. 

The strength of the magnetic field is controlled 
by a unique adjustable parameter Bq that defines 
the absolute value of the loop top magnetic field 
vector B , which is perpendicular on the loop-top 
cross section. Since the vector Bo uniquely deter- 
mines the dipole moment /x = — B i? 3 , the mag- 
netic field vector in any point characterized by the 
position vector r = {x, y, z} relative to the dipole 
origin is given by 

B= 3(Mr)r-rV (1JJ) 
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